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Abstract 

The effect of narrow dibaryon resonances to nuclear matter and structure of 
neutron stars is investigated in the mean-field theory (MFT) and in the relativis¬ 
tic Hartree approximation (RHA). The existence of massive neutron stars imposes 
constraints to the coupling constants of the oj- and cr-mesons with dibaryons. We 
conclude that the experimental candidates to dibaryons di(1920) and d’(2060) if 
exist form in nuclear matter a Bose condensate stable against compression. This 
proves stability of the ground state for nuclear matter with a Bose condensate of 
the light dibaryons. 


The prospect to observe the long-lived H-particle predicted in 1977 by R. Jaffe |T| 
stimulated considerable activity in the experimental searches of dibaryons. It was pro¬ 
posed to examine the H-particle production in different reactions J2J. The experiments 
H did not give a positive sign for the H-particle, however, the existence of the H-particle 
remains an open question which must eventually be settled by experiment. The non- 
strange dibaryons with exotic quantum numbers, which have a small width due to zero 
coupling to the AUV-channel, are promising candidates for experimental searches |4j]. The 
data on pion double charge exchange (DCE) reactions on nuclei 0 exhibit a peculiar 
energy dependence, which can be interpreted || as evidence for the existence of a narrow 
d’ dibaryon with quantum numbers T = 0, J p = 0“ and the total resonance energy of 
2063 MeV. Recent experiments at TRIUMPF (Vancouver) and CELSIUS (Uppsala) seem 
to support the existence of the d’ dibaryon jTJ. A method for searching narrow, exotic 
dibaryon resonances in the double proton-proton bremsstrahlung reaction is discussed in 
Ref. [||. Recently, some indications for a di(1920) dibaryon in this reaction have been 
found [0. 

When density of nuclear matter is increased beyond a critical value, production of 
dibaryons becomes energetically favorable. Dibaryons are Bose particles, so they condense 
in the ground state and form a Bose condensate fL(J |ll[. An exactly solvable model for 


a one-dimensional Fermi-system of fermions interacting through a potential leading to 
a resonance in the two-fermion channel is analyzed in Ref. [12], The behavior of the 


system with increasing the density can be interpreted in terms of a Bose condensation of 
two-fermion resonances. The effect of narrow dibaryon resonances on nuclear matter in 
the mean held theory (MFT) is analyzed in Refs. |[T3| . [TI|| . In the limit of vanishing decay 
width, a dibaryon can be approximately described as an elementary held. 

Despite the dibaryon Bose condensate does not exist in ordinary nuclei, dibaryons 
affect properties of nuclear matter and the ordinary nuclei through a Casimir effect. 
Presence of the background cr-meson mean held inside of nuclei modifies the nucleon and 
dibaryon masses and in turn modifies the zero-point vacuum fluctuations of the nucleon 
and dibaryon helds. This effect contributes to the energy density and pressure. It can 
be evaluated within the relativistic Hartree approximation (RHA). For nucleons, this 
effect is well known |U|. In the loop expansion of quantum hadrodynamics (QHD), MFT 


corresponds to the lowest approximation (no loops), while RHA corresponds to the one- 
loop approximation in a calculation of the equation of state for nuclear matter. 

At zero temperature, a uniformly distributed system of bosons with attractive poten- 

H 


tial is energetically unstable against compression and collapses 


In such a case, the 


long wave excitations (sound in the medium) have imaginary dispersion law: The square 
of the sound velocity is negative a 2 < 0. The amplitude of these excitations increases with 
the time, providing instability of the system. It is necessary to analyze dispersion laws 
of other elementary excitations also. We shall see, however, that in MFT and RHA only 
sound waves can generate an instability. The ground state of nuclear matter with a Bose 
condensate of dibaryons is stable or unstable against small perturbations according as 
the repulsive cu-meson exchange or the attractive a-meson exchange is dominant between 
dibaryons. 

In this paper, we investigate the hypothesis that the dibaryon matter is unstable 
against compression. In such a case, formation of dibaryons in nuclear matter can be 
treated as a possible mechanism for a phase transition into the quark matter. If central 
density of a massive neutron stars exceeds a critical value for formation of dibaryons, the 




















neutron star should convert into a quark star, a strange star, or a black hole. Some of the 
observed pulsars are identified quite reliably with ordinary neutron stars [lTj]. From the 


requirement that the dibaryon formation is not energetically favorable at densities lower 
than the central density of neutron stars with a mass 1.3M 0 , we derive constraints to 
the coupling constants of the mesons and dibaryons di(1920) and d'(2060) and conclude 
that narrow dibaryons in this mass range can form a Bose condensate stable against 
perturbations only. The effect of the dibaryons to stability and structure of neutron stars 
in different phenomenological models is analyzed in Refs, 
binding energy of strange matter 


11, 181. Constraints to the 


191 from the existence of massive neutron stars are 


discussed in Ref. 20 


The dibaryonic extension of the Walecka model |15|] is obtained by including dibaryons 
to the Lagrangian density 


£ = ^{iduTfn - m N - g a cr - g u uj^)^ + 


I 2 2 


\F 2 + + {dp - + iKu^ip - (m D + h a a) 2 ip*ip. 


( 1 ) 


Here, T is the nucleon field, u; M and a are fields of the u- and cr-mesons, F^ u = d v uj jJ — d ll uj lJ1 
ip is the dibaryon isoscalar-scalar (or isoscalar-pseudoscalar) field. The values and 
m a are the u- and cr-meson masses and the values <£, g a , h w , h„ are coupling constants 
of the u- and cr-mesons with nucleons (g) and dibaryons ( h ). 

The (T-meson mean field a c determines the effective nucleon and dibaryon masses in 
the medium 


^N T 9a@ci 
m* D = m D + h a a c . 

The nucleon scalar density in the RHA is defined by expression |T3 


Pns =< t(0)T(0) >= 7 


dp 


m 


N 
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■0(Pf - |p|) - 4m 3 N t(m* N /m N ) 


( 2 ) 

( 3 ) 

( 4 ) 


where 


47 t 2 ((x) = x 3 hix + 1 


x-\( l -x ) 2 + y^-A) 3 - 


The last term in Eq.(4) occurs after the renormalization of the scalar density. Here, 7 = 2 
for neutron matter and 7 = 4 for nuclear matter. 

We investigate here properties of the nuclear matter below the critical density for 
formation of dibaryons, so the dibaryon condensate is zero < |</?(0)| >= 0. The vacuum 
contribution to the scalar density of dibaryons can be found to be 


2 m* D p DS = 2 m* D < £(0)V(0) >= m 3 D ((m* D /m D ). 


( 5 ) 


It differs from the nucleon term in the sign and in the statistical factor (one should replace 
4 (= 2 S x 2/) —> 1). The self-consistency condition for the nucleon effective mass has the 
form 

m* N = m N - ^(gapNS + h s 2m* D p DS ). (6) 

Tflrr 


The renormalized vacuum contribution to the nucleon energy-momentum tensor has 
the form [|15|] 

~*N ((\\ \ _ a ^4 ^^ \ 


< >vac= -4g^rn%rj(m* N /m N ) 






































where 


7 IS or. 

167r 2 p(x) = x 4 /nx + 1 — x — -(1 — x) 2 + —(1 — x) 3 — —(1 — x) 4 . 

Z O _L Zj 

For dibaryons, we get expression 

< T%,(0) >vac= g^m A D ri(m* D /m D ). 


( 8 ) 


The elementary excitations in nuclear matter with a Bose condensate of dibaryons 
correspond to nucleons and antinucleons, cn-mesons, cr-mesons, and dibaryons and an- 
tidibaryons. The dispersion laws for these quasiparticles are found in Ref. JOj]. The 


nucleon and antinucleon dispersion laws have the same form as in the vacuum with a 
replacement tun —> m* N and therefore cannot generate an instability of the system. The 
dispersion laws of cu-mesons, cr-mesons, and antidibaryons turn out to be real also. The 
possible source of the instability are the dibaryon quasiparticle excitations only, which 
are responsible for a long wave perturbations of the system and connected with existence 
of the sound in the medium. 

The square of the sound velocity has the form [|hj 


where 


9 a 
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We see that a 2 s > 0 for 


h 2 h 2 

n uj n a 

ml 77il' 

(11) 


The validity of inequality (11) is sufficient condition for stability of the ground state of 
nuclear matter with a Bose condensate of dibaryons. 

The physical meaning of the inequality (11) can be clarified by considering the inter¬ 
action energy of a uniformly distributed dibaryon matter pw( x i) = Pdv = constant: 


W = ^ j dxidx 2 pw(xi)pw(x 2 )R(|xi - x 2 |). 

The Yukawa potential V (r) for two dibaryons has the form 

' hi e- m ” r 
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The integration gives 
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(14) 


where Nd is the number of dibaryons. When the dibaryon density p^y increases, the 
energy for a 2 > 0 increases also, the pressure is positive, and so the system is stable. 

At present, the coupling constants of the mesons with dibaryons are not known 
with precision good enough to draw a definite conclusion concerning the stability of the 
dibaryon matter. 






Here, we show that violation of the inequality (11) for light dibaryons is in contradic¬ 
tion with the existence of massive neutron stars. 

Given that the ratio h a / (2 g a ) between the cr-meson couplings with dibaryons and 
nucleons is fixed, one can find the c o- and cr-meson couplings with nucleons by fitting the 
nuclear matter binding energy E/A — m n = —15.75 MeV at the empirical equilibrium 
density po = 0.148 fm -3 determined from the density in the interior of 20S Pb |H|. It 
corresponds to the equilibrium Fermi wavenumber kp = 1.3 fm -1 . The dependence of 
the effective nucleon mass at the equilibrium density on the ratio h a / (2 g a ) is shown 
on Fig.l (a). In Figs.l (b) and (c) we give dependence of the incompressibility K = 
9p 0 (d 2 £/dp 2 )\ p=po and the asymmetry coefficient a 4 on the ratio h a / (2 g a ). In Fig.l (d), the 
values C 2 S = g 2 (m,N/vn a ') 2 and C 2 = g 2 (rripr jrn w ) 2 are plotted. Notice that h a /(2g a ) = 0 is 
equivalent to RHA with no dibaryons. For comparison, we give the results of MFT where 
the effect of dibaryons to the nuclear matter below the critical density for occurrence of a 
Bose condensate of dibaryons is absent. The results RHA for different dibaryons are not 
much different. In this work, we study the physical implications of the effective Lagrangian 
(1) describing nucleon and dibaryon degrees of freedom. Other baryons can be included 
in the QHD framework in a similar way. Their effect is quite small. We have checked that 
the Casimir effect caused by the inclusion of other octet baryons (2 x 8 = 16 degrees of 
freedom) shift the critical value h„ only by about 25% if one assumes that the cr-meson 
is an SU(3)/ singlet. The inclusion of octet and decuplet baryons (2 x 8 + 4 x 10 = 56 
degrees of freedom) with a universal sigma-meson coupling constant increases the critical 
value of h a /(2g a ) to 1.2. 

When the ratio h a / (2 g a ) approaches a value 0.8, the system of equations blows up and 
the empirical equilibrium properties of the nuclear matter can no longer be reproduced. 
When x —> 0 ((x) = 0(x 4 ), so the zero-point contributions to the scalar density of 
nucleons and dibaryons, which have the opposite signs, are comparable for Ag 4 /rriN ~ 
/ij/rai). The dibaryon effects become large for h a /(2g a ) & 0.5(A?np,/mjv) 1 ^ 4 ~ 0.84. The 
greater dibaryon mass, the greater the upper limit to the ratio h a /{2g a ). This effect is 
seen in Fig.l. 

The saturation curve is shown in Fig.2 and the effective nucleon mass dependence 
on the Fermi wavenumber kp is shown in Fig.3 for h cr /(2g 0 -) = 0.6 in case of the H- 
particle. EOS in RHA is softer than in MFT. The contributions of the vacuum zero-point 
fluctuations of nucleons and dibaryons partially cancel each other, so the inclusion of the 
dibaryons makes the EOS stiffer. In Fig.2, we see that the dashed curves corresponding to 
RHA with dibaryons lie above the solid lined corresponding to the RHA with no dibaryons. 
The same effect is seen on Fig.3. In MFT, the nucleon effective mass decreases with the 
density faster then in RHA. Due to the partial compensation of the nucleon and dibaryon 
contributions to the vacuum scalar density, the dashed lines lie below the solid lines. 

In Fig.4 we show the critical densities for occurrence of dibaryons H(2220), d’(2060), 
and di(1920) in MFT and RHA in nuclear and neutron matter versus the ratio h a / (2 g a ) for 
huj/h™ ax = 1, 0.8, and 0.6 where h™ ax = h a m u lm n is the maximum value for the cu-meson 
coupling constant with dibaryons at which the inequality (11) is violated (h u /h™ ax = 1 
corresponds to a 2 = 0 and h u /h™ ax = 0.8 and 0.6 correspond to a 2 < 0). In RHA, 
dibaryons occur at higher densities. The coupling constant h u determines the energy 
of dibaryons in the positive uxmeson mean field. The greater the h^, the greater the 
density is required to make production of dibaryons energetically favorable. This effect 
is seen in Fig.4: The solid lines h UJ /h™ ax = 1 lie above the long-dashed and dashed lines 



hu,/h™ ax = 0.8 and 0.6, respectively. 

When the dibaryon matter is unstable against compression, production of dibaryons 
with increasing the density results to instability of neutron stars with subsequent phase 
transition into the quark matter and conversion of neutron stars into quark stars, strange 
stars, or black holes. In such a case, the maximum masses of neutron stars are determined 
by the mass and the coupling constants of the mesons with the lightest dibaryon. In Fig.5 
we show the minimal neutron star masses in which dibaryons can occur. 

The MFT and RHA EOS for neutron matter at supranuclear densities are matched 
smoothly with the BBP EOS f2]| at densities pd r i P < p < 0.8 p 0 where pd r i P = 4.3 
10 n g/cm 3 and then with the BPS EOS [22| at densities p < Pdri P ■ The maximum neu¬ 


tron star masses are sensitive to the value of the equilibrium Fermi wavenumber. If we 
chose kp = 1.42 fm -1 instead of kp = 1.3 fm -1 , the maximum masses in MFT (with no 
dibaryons) are reduced from 3 M & down to 2.6M 0 [|H|. The choice kp = 1.3 fm -1 provides 
less stringent (therefore more conservative) constraints to the meson-dibaryon coupling 
constants. We do not show the results for the di(1920) dibaryon, since its condensation 
starts at a density p ~ Po providing conversion to quark stars of neutron stars with very 
low masses M < 0.2 M & . 

In Fig. 6 we show the parameter space for the coupling constants of dibaryons with 
the mesons. Our discussion and the validity of our arguments are restricted to the region 
a 2 < 0 in which the dibaryon matter is unstable against compression. The requirement of 
stability of the normal nuclear matter at the saturation density allows to get constraints to 
the coupling constants. The corresponding curves (straight lines in the MFT case) marked 
by arrows with the white ends restrict from below the parameter space of the coupling 
constants. The dotted line in the MFT case, which refers to the di(1920) dibaryon, is 
very close to the dashed-dotted line a 2 = 0. For low values h a the dotted line lies above 
the line a 2 = 0. In such a case, the dibaryon matter unstable against compression cannot 
exist. The window in the parameter space for the unstable dibaryon matter for di(1920) 
is, however, much greater for RHA. 

If a conservative assumption is used, namely, that pulsars with a mass 1.3M 0 are 
ordinary neutron stars, the constraints to the meson coupling constants with dibaryons 
can further be improved. The corresponding curves (straight lines in the MFT case) are 
shown on Fig.6. We see that the dotted and dashed lines lie above or very close to the 
line a 2 = 0. It means that for di(1920) and d’(2060) dibaryons, the dibaryon matter 
unstable against compression cannot exist (owing to a very small window for d’(2060) at 
higher values of the h a ). For the H-particle, there is a window in the parameter space 
between the line a 2 = 0 and the solid curves marked by arrows with the black ends, which 
corresponds to the dibaryon matter unstable against compression. 

The H-particle interaction were studied in the non-relativistic quark cluster model 


23, 24] successful in describing the AW-phase shifts. The coupling constants of the mesons 


with the H-particle can be fixed by fitting the depth and the position of the minimum 
of the HH-adiabatic potential ||25| to give h u /(2g 0J ) = 0.89 and h a /(2g a ) = 0.80. These 
values are marked on Fig.6 with a cross. These estimates are used in the MFT-calculations 
[T3[ |T4|j . They correspond to the unstable dibaryon matter and are in the allowed region 
of the parameter space for the H-particle. The energetically favorable compression of the 
H-matter can lead to the formation of the absolutely stable strange mater |T9[ , producing 


conversion of neutron stars to strange stars 0 


MFT and RHA EOS both are very stiff. In the soft models like the Reid one (for 































a review of the nuclear matter models see 0), the central density of neutron stars is 
much greater than in stiff models. Respectively, conditions for occurrence of new forms of 
nuclear matter are more favorable. From Figs. 5 and 6, we see that the softer RHA EOS 
produces lower upper limits to the neutron star masses and, respectively, more stringent 
constraints to the meson-dibaryon coupling constants as compared to the stiffer MFT 
EOS, despite in RHA dibaryons occur at higher densities (see Fig.4). One can assume 
that this effect is of the general validity and that softer EOS like the Reid one give more 
stringent constraints to the meson-dibaryon coupling constants. We consider therefore 
the constraints given in Fig.6 as the conservative ones. 

In the conclusion, we showed that the hypothesis of instability of the dibaryon matter 
against compression is in contradiction with the hypothesis that pulsars of a mass 1.3M 0 
are ordinary neutron stars for dibaryons di(1920) and d’(2060). This conclusion is valid 
for all narrow dibaryons with the same quantum numbers in the same mass range. The 
H-particle is sufficiently heavy, its condensation starts at higher densities, respectively, 
constraints to the meson-dibaryon coupling constants are not much stringent, allowing 
a possibility for the H-particle to form a condensate unstable against compression with 
subsequent formation of the strange matter. The meson coupling constants with dibaryons 
di(1920) and d’(2060) should obey the inequality (11). The meson coupling constants with 
the H-particle should lie above the solid curves on Fig.6. 
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Figure captions 


Fig.l The effective nucleon mass (a), the asymmetry coefficient (b), the incom¬ 
pressibility of nuclear matter (c) at the saturation density, and the coupling constants 
C'l = g 2 s {rriN / m s ) 2 and C 2 = g^rn^ / m,J) 2 in the RHA versus the ratio h rJ /(2g rJ ) between 
the a-meson coupling constants with dibaryons and nucleons. The MFT results are shown 
for comparison. The coupling constants are fixed by fitting the minimum and depth of the 
energy per baryon number at the saturation density of nuclear matter. The solid, dashed, 
and dotted curves correspond to the dibaryons H(2220), d , (2060), and di(1920). The 
reasonable description of the properties of the nuclear matter at the saturation density is 
possible for h a /(2g a ) < 0.8. 

Fig.2 Saturation curve for nuclear matter in RHA with no dibaryons (solid line) and 
in RHA with inclusion of the H-dibaryon (dashed line) for h a /(2g a ) = 0.6. 

Fig.3 The effective nucleon mass versus the fermi momentum of nucleons in the nuclear 
matter (upper curves) and in the neutron matter (lower curves). The solid lines correspond 
to RHA with no dibaryons, the dashed lines correspond to RHA with inclusion of the H- 
dibaryon. 

Fig.4 Critical densities for occurrence of the dibaryons H(2220), d’(2060), and di(1920) 
in nuclear and neutron matter in MFT and RHA versus the cr-meson coupling constant 
with dibaryons for hu/h™ ax = 1, 0.8, and 0.6 (the solid, longed-dashed, and dashed lines, 
respectively), where h™ ax = h a m^lm a is the maximum value for the uj- meson coupling 
constant with dibaryons at which the dibaryon matter is unstable against compression 
(see the text). 


Fig.5 The minimum neutron star masses, in which the dibaryon formation becomes en¬ 
ergetically favorable, versus the a-meson coupling constant with dibaryons for h^/h^ ax = 
1 , 0.8, and 0.6 (the solid, longed-dashed, and dashed lines, respectively) where h™ ax is 
the maximum value for the cu-meson coupling constant at which the dibaryon matter is 
unstable against compression (a 2 < 0). The results are given in MFT and RHA and for 
two dibaryons H(2220) and d’(2060). The minimum neutron star masses for the di(1920) 
are very small (< 0.2M 0 ). 


Fig.6 Parameter space for the coupling constants of the a- and uj- mesons with dibaryons 
in MFT and RHA. The dashed-dotted line a 2 = 0 divides the parameter space into two 
parts. The upper left part of the parameter space corresponds to the dibaryon matter 
stable against compression (square of the sound velocity is positive a 2 > 0), the lower 
right part corresponds to the nuclear matter unstable against compression (a 2 > 0). The 
solid, dashed, and doted curves constrain from below the regions in which the dibaryon 
formation is energetically not favorable in ordinary nuclei (curves marked by arrows with 
white ends) and at the density equal to the central density of a 1.3M 0 mass neutron star 
(curves marked by arrows with black ends). The cross refers to the H-particle coupling 
constants with the mesons, determined from the adiabatic potential |25| . 



